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I. INTRODUCTION 

The simulation of rheological properties of molecular 
fluids requires realistic temperature control correctly tak- 
ing into account the streaming velocity of the fluid. This 
has become feasible following the introduction of ther- 
mostatting techniques defining temperature using only 
configurational variables. Configurational approaches are 
also important in many other contexts, for example in 
systems subject to an external pressure constraint and in 
understanding conformational dynamics of biomolecules. 

In a series of papers^— Braga and Travis proposed a 
Nose-Hoover type thermostat based on a configurational 
definition of temperature (in short BT thermostat), re- 
moving difficulties associated with the earlier Delhomelle- 
Evans approach^ which for example does not preserve the 
canonical distribution in the equilibrium case. Recently 
the BT thermostat has had further attention, develop- 
ment and application in the literature^—. In the prin- 
cipal paper—, the authors considered the configurational 
thermostat scheme proposed in reference o 12 i 13 (in short 
SDC scheme) as a starting point and then presented their 
own configurational temperature thermostat as a further 
development of the SDC scheme. Correct action of both 
BT and SDC thermostats has been confirmed by a vari- 
ety of numerical simulations in these works. 

In this note, we link together two above-mentioned 
contributions in the development of configurational ther- 
mostats and thus accomplish two goals. First, we show 
that the configurational temperature BT thermostat is a 
particular case of the SDC scheme. Second, we demon- 
strate how this particular case benefits from the general 
metho d 12 ' 13 in a variety of aspects including the physi- 
cal context, a stochastic counterpart thermostat (to the 
authors' knowledge this has not previously appeared in 
the literature) and ergodicity issues. 



II. THE CONFIGURATIONAL TEMPERATURE 
THERMOSTAT 

Let us clearly state in what way the BT thermostat is 
a particular case of the SDC scheme. The general SDC 
schem o 12 ' 13 involves two different definitions of the con- 
figurational temperature and the corresponding double 



temperature control (for example, to allow control at two 
timescales), coupling between thermostats, anisotropy, 
and a series of definitions of the thermostat variable that 
does not directly relate to a particular temperature con- 
trol. To reduce this general configurational thermostat 
method to the BT equations (our first goal) , consider the 
simplified version of the SDC where the virial tempera- 
ture control and the coupling of thermostat variables are 
omitted: 
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Then consider the dynamic variables £k according to the 
case (a) in reference d 12 ' 13 (variables £ for all particles of 
the system are varied independently), 

ik = --^-—v^viq), 

set = 1, and finally define the canonical variables, 
r fc = \/rak(\k, Pfc = \/m>k£,k- As a result, we arrive at 
the BT equations as a particular case of the more general 
SDC scheme, 
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easy to show that the canonical distribution, 

exp[— - k -^p{^2 2m V)]> 1S invariant for this dynamics. 

In the following paragraphs we demonstrate in what 
way the SDC scheme contributes to the understanding 
of the configurational temperature thermostats. 

A frequent problem, even if the canonical measure is 
invariant for the dynamics ([T]), is that only ergodic dy- 
namics correctly samples the canonical distribution. In 
general, this condition is very difficult to prove for sys- 
tems with many degrees of freedom. In particular, simple 
one-dimensional systems^ reveal a difficulty and provide 
an important challenge for deterministic thermostatting 
methods. On the other hand, stochastic dynamics is typ- 
ically ergodi o 15 ' 16 . 



III. THE STOCHASTIC COUNTERPART 
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The Nose-Hoover thermostat is a deterministic coun- 
terpart of the Langevin stochastic dynamics (e.gr^). 
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This fact is useful to determine the appropriate charac- 
teristic time scale of the dynamics. To clarify the physical 
sense of the configurational temperature thermostat, we 
formulate a stochastic counterpart of the dynamics ([1}. 
For the sake of brevity we provide only the stochastic 
dynamics and its justification. A formal analysis can be 
done following the method presented in reference a 17 ' 18 . 

A particle of mass m in a dissipative medium under 
the influence of a constant force F, has achieved the con- 
stant velocity v s — tF, where r is the mobility, after the 
transient time, ~ mr. We take into account the mobility 
as a perturbation of velocities. Explicitly, consider the 
following stochastic equations of motion: 

q fe = ^ - rV qfc y(g) + V2Dh(t), 
m k 

Pk = -V q ^(<z), (2) 

where {fk(t)} is the set of independent standard "white 
noises". The Fokker-Planck equation corresponding to 
@ has the form, 

dtP = Y / {-^^[(—-rV (lk V)p} 

(fe) mk 

+ V Pfc • (V^Vp) + DA^p} ee F*p, (3) 

where F* is the Fokker-Planck operator. 

These stochastic dynamics can be generalized to the 
case of the state dependent mobility, r — > T(p(p), where 
tp(p) is a positive function of momenta. The Fokker- 
Planck equation dtp — F*p corresponding to that case 
has the form ([3]) where r — > np(p) and D — > D<p(p). 
Indeed, suppose r and D are connected to each other 
by the relation D = rksT. Substituting the canoni- 
cal distribution, p ~ exp[-/3(^(fe)(2TOfc) _1 p 2 + V (<?))], 
into F*p and then applying the condition, D = rksT, 
we arrive at the identity, F*p = 0. Thus the canonical 
measure, dp ~ Pll(fc) dq k dpk is invariant for this dy- 
namics and we obtain the stochastic counterpart of the 
deterministic configurational thermostat (|TJ) . The stan- 
dard case is tp(p) = 1, F*p = 0. The more general case, 
<p = tp(p, q), can be considered but in that case we must 
specify the interpretation of the stochastic equations (Ito 
or Stratonovich) and correspondingly add an extra drift 
term into {2]). 

IV. THE EXTENDED CONFIGURATIONAL 
THERMOSTAT 

Following the general SDC schem e 12 ' 13 the determin- 
istic BT thermostat ((TJ) can be refined in the following 
manner. Consider the extended dynamical system: 

h = — - T~V Tk V + £fc, p k = -Vr k V 1 

1 N 

t = Y) — [(Vr k V) 2 - k B TA rk V] , £ k - h fe , (4) 
^ TT k=i 



where {£k} are additional dynamical variables. Func- 
tions {hj.} can be fixed in a variety of forms as described 
in reference o 12 ' 13 . For example, set hfe = —-^J2(i) ^r { V 
for all particles of the system. Then the canonical mea- 
sure (augmented with the Gaussian measure for thermo- 
stat variables r and £) is invariant for dynamics (|4]). In 
that case, the velocity of the center of inertia of a system 
dynamically fluctuates around zero with mobility Q^ 1 , 
where is a parameter. 

V. ON ERGODICITY 

The thermostat dynamics is able to correctly sample 
the canonical distribution only if the motion is ergodic. 
To improve ergodicity, we proposed in reference^ 3 - to per- 
turb the dynamics of the augmented thermostat vari- 
ables (the Gaussian invariant measure for these variables 
is supposed) by additional random noise. This method 
was further investigated together with mathematical jus- 
tification in reference o 15 ' 16 . Numerical effectiveness of 
this method was also confirmed— ' 1 5 1 1 6 . Taking into ac- 
count additional advantages of this method including the 
"gentle" action of the stochastic ter m 15 ' 16 , we propose to 
apply it to configurational temperature thermostats ([1} , 
(|4|) similarly to the description in reference^ 3 -. For ex- 
ample, in set ([1]) we propose to perturb the dynamics of 
t- variable by the "white noise" f(t), 

i N 

f =7rT, i^r k V) 2 - k B TA rk V] - At + s/ZDf(t), 
Wtt fe=i 

where (positive) A and D are connected to each other by 
the relation, k B TA — DQ TT . Then the augmented canon- 
ical measure is invariant for the dynamics. We conjec- 
ture that this is the unique invariant measure (compare 
withi 5 -). 

VI. CONCLUSION 

In this Note, we have linked together two contribu- 
tions in the development of configurational thermostats, 
the BT thermostat and the SDC scheme. We have shown 
that proposed configurational thermostatsir— are gener- 
alized and enriched in both understanding and content 
by the SDC scheme. We have presented the stochastic 
counterpart ((2|) to the configurational thermostat ([TJ. 
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